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Abstract 

Using auxiliary-mass method, O(N) invariant scalar model is in- 
vestigated at finite temperature. This mass and an evolution equa- 
tion allow us to calculate an effective potential without an infrared 
divergence. Second order phase transition is indicated by the effective 
potential. The critical exponents are determined numerically. 
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Symmetry restoration of 0(N) scalar model at high temperature is very 
important since many physical systems belong to the same universality class 
as it: polymer phase transition (N— ►()), critical liquid- vapor phase transi- 
tion (N=l), alloy (e.g. /3-brass) phase transition (N=l), uniaxial ferromag- 
net phase transition (N=l), superfluid phase transition (N=2), ferromagnet 
phase transition (N=3), and chiral phase transition with two flavor massless 
quarks (N=4) 0, g. 



The phase transition had better be investigated by finite-temperature 
field theory, which is based on the statistical principle only. The perturba- 
tion theory, however, breaks down around the critical temperature when the 
phase transition is second or weakly first order ||, even if daisy- diagrams 
(ring-diagrams) [§],§, |6[ are resumed. Investigation into phase transitions at 
finite-temperature has long been hampered by this failure. Many methods to 
avoid the failure were proposed: CJT method [[!]], renormalization improve- 
ment H, novel summation ||, pade improvement fllQfl , exact renormalization 
group with temperature 



TTJ , and auxiliary-mass method [T2|, [13], 



The auxiliary-mass method is used in the present paper. The method is 
based on the following idea. First an effective potential can be calculated at 
large mass by the perturbation theory since it is reliable there. Second the 
effective potential is extrapolated to the small mass range, where the per- 
turbation theory is not reliable, using an evolution equation. Finally various 
quantities (e.g. critical exponents) are determined from the effective poten- 
tial. The method is applied to O(N) scalar model concretely in the following; 
the phase transition of the model is investigated and the critical exponents 
are determined. 



We consider the following Lagrangian density, 

£* = -\{jt)' ' ~\^^f -^1-^1? + ^ + c± , (1) 

where J a is an external source function and the index "a" runs from 1 to 
N. First the effective potential V is calculated by the perturbation theory 
within one- loop order at large mass m 2 = M 2 = 0(T 2 ). We choose field 
expectation values (fi a = (p 5\ a without a loss of generality because of O(N) 
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We note that the daisy-resummation is not necessary because of the large 
mass and one-loop zero-temperature effect is negligible if the coupling is 
weak. Next this effective potential is extrapolated to smaller mass using the 
following evolution equation [O], 
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Here IT^ = U Sig (m 2 , <p 2 , p 2 , p 2 , T) and U NG = U NG (m 2 , 2 , p 2 , p 2 , T) are the 
full self-energy, which are those for massive and massless mode in the broken 
phase. This equation is modified from that of A0 4 theory [[13] straitforwardly 
through a diagonalization of the propagator. Though this equation is exact, 
it can not be solved without an approximation; because it includes the full 
propagator which can not be known exactly. We then replace as follows, 
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This replacement corresponds to fixing an external momentum of the full 
self-energy to zero 1 . It allows us to convert eq.(§) to the following partial 
differential equation, 
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1 This is the first approximation of a systematic calculation 15 
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Figure 1: The effective potential obtained by the auxiliary-mass method 
(N = 4, A = 1). Second-order phase transition occurs at the critical temper- 
ature. Similar behavior are observed for the other N and A. 
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We numerically solve the equation to m 2 = —fi 2 , where O(N) symmetry is 
broken at zero-temperature, under the initial condition eq.(^). We solve it 
above the critical temperature only and can get sufficient information about 
the critical phenomenon. 



The effective potential are calculated using numerical methods in Ref . [|I~4 



They are shown in fig.JTJ. We can not observe a second minimum which is 
a hallmark of a first order phase transition. Instead, the curvature of the 
origin decreases and vanishes smoothly as temperature decreases; one can, 
therefore, observe that the phase transition of the model is second order as 
it should be. 



The critical temperature is shown as a function of N in fig.0. Though 
they resemble to a result of the perturbation theory at leading order 
T, 
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which is determined from the condition that mass with the 



A(Af+2) 

daisy-diagram vanishes, they have a slight difference quantitatively. 
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Figure 2: Critical temperature as a function of N at A = 1 (o). This resembles 

to a result of the perturbation theory at leading order ( ) but has a slight 

difference quantitatively. 



Finally the critical exponents, 7 and 5, are determined using this poten- 
tial. We determine 7 from a second derivative of the effective potential with 
respect to at the origin. The critical exponent 7 is defined as follows, 



X 



dJi 



[t=(T-T c )/T c ) 



(6) 



Ja=0 



the following identity then relates the second derivative to the susceptibility, 
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Here <p c is determined from the condition, ~ _ =0. The second derivative 

0<t> <j>=(j> c 

are shown in log-scale in fig.|3]. We then determine the gradients, which is the 
very 7 we want. One can observe that they become steeper as N increases; 
then, 7 become larger. The results are 7 = 1.13 (N — > 0), 1.37 (N = 
1), 1.47 (N = 2), 1.60 (N = 3), 1.66 (N = 4). These results are summarised 
in table. ([I]), compared with the Landau approximation and world best values. 
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Figure 3: Second derivative of the effective potential with respect to (p. The 
gradients are steeper for larger N. 



We determine S using the effective potential at the critical temperature. 
The critical exponent 5 is defined as follows. 
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The following relation, which derived from eq.([8]), enables us to determine S 
from the effective potential at the critical temperature, 
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The effective potential at the critical temperature is shown in fig.[3]. We de- 
termine 5 from it. The result are 5 = 3.8 (N — > 0), 4.0 (N = 1), 4.2 (N = 
2), 4.4 (N = 3), 4.4 (N = 4). 



In the present paper, we showed that the phase transition of O(N) scalar 
model is second order using auxiliary-mass method at finite temperature. 
This is a great progress because we can not show it using the perturbation 
theory with daisy-resummation - - the traditional method to calculate the 
effective potential at finite temperature [[|. 

Since the phase transition turned out to be second order, we determined 
the critical exponents of the phase transition. Though the results are not as 
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7 (LA,WBV) 


5 (LA,WBV) 




1.13 (1, 1.16) 


3.8 (3, 4.77) 


N=l | 
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1.37 (1, 1.24) 


4.0 (3, 4.76) 


N=2 | 


1 


1.47 (1, 1.32) 


4.2 (3, 4.78) 


N=3 | 


1 


1.60 (1, 1.40) 


4.4 (3, 4.78) 


N=4 | 




1.66 (1, 1.48) 


4.4 (3, 4.85) 



Table 1: The critical exponents, 7 and 5, obtained in the present paper. 
Those of Landau approximation (LA) and world best values (WBV) are also 
summarised. We used lattice results as WBV. 

accurate as world best values, they are much better than that obtained in 
the Landau approximation. The error would be due to the replacement (f|). 
We must improve the approximation in order to get more accurate values [|T5|]. 

In conclusion, the auxiliary-mass method turns out to be reliable not 
only qualitatively but also in rough quantitative estimation in O(N) invari- 
ant scalar model. This method would be reliable in the other models. What 
is more, we can investigate not only second-order phase transition but also 



first-order phase transition without any modifications [ I7| . This method, 
therefore, enable us to investigate the phase transition of various models: the 
cubic anisotropy model, the abelian Higgs model, and the standard model. 

We finally express our thanks to T. Inagaki for valuable discussions and 
communications. J.S. is supported by JSPS fellowship. 
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